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We study the nonequilibrium time evolution of the Bose-Einstein condensate of interacting bosons 
confined in a leaky box, when its number fluctuation is initially [t — 0) suppressed. We take account 
of quantum fluctuations of all modes, including k = 0, of the bosons. As the wavefunction of the 
ground state that has a definite number TV of interacting bosons, we use a variational form \N,y}, 
which is obtained by operating a unitary operator e*''^^' on the number state of free bosons. Using 
e'^'fy), we identify a "natural coordinate" So of the interacting bosons, by which many physical 
properties can be simply described. The |A'^, y) can be represented simply as a number state of 6o, 
we thus call it the "number state of interacting bosons" (NSIB). To simulate real systems, for which 
if one fixes A*' at t = will fluctuate at later times because of a finite probability of exchanging 
bosons between the box and the environment, we evaluate the time evolution of the reduced density 
operator p{t) of the bosons in the box as a function of the leakage flux J. We concentrate on 
the most interesting and nontrivial time stage, i.e., the early time stage for which Jt <C A, much 
earlier than the time when the system approaches the equilibrium state. It is shown that the time 
evolution can be described very simply as the evolution from a single NSIB at t < 0, into a classical 
mixture, with a time dependent distribution, of NSIBs of various values of A' at t > 0. Using bo, 
we successfully define the cosine and sine operators for interacting many bosons, by which we can 
analyze the phase fluctuation in a fully quantum-mechanical manner. We define a new state N, y) 
called the "number-phase squeezed state of interacting bosons" (NPIB), which is characterized by 
a complex parameter ^. It is shown that p{t) for t > can be rewritten as the phase-randomized 
mixture (PRM) of NPIBs. Among many possible representations of p{t), this representation is 
particularly convenient for analyzing the phase fluctuations and the order parameter. We study the 
order parameter according to a few typical definitions, as well as their time evolution. It is shown 
that the off-diagonal long-range order (ODLRO) does not distinguish the NSIB and NPIB. Hence, 
the order parameter H defined from ODLRO, does not distinguish them, either. On the other hand, 
the other order parameter ^, defined as the expectation value of the boson operator ^p, has different 
values among these states. In particular, for each element of the PRM of NPIBs, we show that 
^l/ evolves from zero to a finite value very quickly. Namely, after the leakage of only two or three 
bosons, each element acquires a full, stable and definite (non-fluctuating) value of 

PACS numbers: 03.75.Fi, 05.30.Jp, 05.70.Fh, 05.70.Ln 



I. INTRODUCTION AND SUMMARY 

The Bosc-Einstcin condensation (BEC) has been observed in various systems 0], including the liquid Helium 
excitons in photo excited semiconductors , and atoms trapped by laser beams . Although BEC was originally 
discussed for free bosons, a condensate of free bosons does not have the superfluidity pi, hence many-body interactions 
are essential to interesting behaviors of condensates. The condensed state of interacting bosons in a box of finite volume 
V is conventionally taken as the state in the Bogoliubov approximation, which we denote by |ao, y'^')'^' [see Refs. 
and Eq. ( |2C| ) below]. In this state, the boson number N has finite fluctuation, whose magnitude is {6N'^) ^ (N) [Eq. 
(p9|)]. This fluctuation is non- negligible in small systems, such as Helium atoms in a micro bubble |^ and laser-trapped 
atoms where (N) is typically 10^ - 10'', and thus ^/{SWj/{N) = 3 - 0.1%. This means that in such systems 

if one fixes N with the accuracy better than 3 — 0.1%, then {SN'^) < {N), thus the state |ao,y''')''' is forbidden, and 
another state should be realized. In most real systems, there is a finite probability of exchanging bosons between the 
box and the environment. Hence, if one fixes N at some time, N will fiuctuate at later times. Namely, the boson 
state undergoes a nonequilibrium time evolution when its number fluctuation is initially suppressed. The purpose of 
this paper is to investigate the time evolution of the interacting bosons in such a case, and to discuss how an order 
parameter is developed. 

We first review and discuss the case where the box is closed and the boson number N is exactly fixed (section O) . 
The ground-state wavefunction of such a case may be obtained by the superposition of Bogoliubov's solution |ao, y/'^^ 
over various values of the phase of ao [Eq. (|3^)]. The resulting state |A^, y'=')'=' has the same energy as jao, y'^')'^' because 
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of the degeneracy with respect to the phase of aQ. (This degeneracy leads to the symmetry breaking.) However, such 
an expression is not convenient for the analysis of physical properties. To find the ground state in the form that is 
convenient for analyzing physical properties, we derive an effective Hamiltonian H [Eq. (||)], which includes quantum 
fluctuations of all modes including k = 0, from the full Hamiltonian Hb of interacting bosons. Here, we neglect effects 
due to spatial inhomogeneity of the boson states in the box, because we are not interested in such effects here, and 
also because we expect that main physics of the nonequilibrium evolution of our interest would not be affected by 
such effects. A renormalization constant Zg appears in H. Although Zg formally diverges [Eq. (Q)], the divergence is 
successfully renormalized, i.e., the final results are independent of Zg and finite. As the ground state of 7?, we use a 
variational form |iV, y), which is similar to that of Girardeau and Arnowitt This form takes a compact form [Eq. 
(^)]: it is obtained by operating a simple unitary operator e**^'^''^ on the iV-particle state of free bosons, where G{y) 
is a simple bi-quadratic function of the bare operators [Eq. (^^]. This state has the same energy as |ao,y^')^', or, 
equivalently, | A^, y"^')"^'. (Precisely speaking, they have the same energy density in the macroscopic limit, i.e., when 
V —^ oo while keeping the density n finite.) 

Using the unitary operator e*'^'^''^ we then identify a "natural coordinate" 60 [Eq. (^7|)] of the interacting bosons, 
by which many physical properties can be simply described (section HI). Unlike the quasi-particle operators obtained 
by the Bogoliubov transformation, 60 is a nonlinear function of bare operators. Moreover, the Hamiltonian is not 
diagonal with respect to bo. Such a nonlinear operator, however, describes the physics very simply. For example, 
\N,y) is simply represented as a number state of bo. We thus call |A, y) the "number state of interacting bosons" 
(NSIB). We can also define, through 60, the cosine and sine operators for interacting bosons [see below]. Moreover, 
using bo, we decompose the boson field ip into two parts [Eq. (|5^)]: one behaves anomalously as ^ 00 and the 
other is a normal part. In the decomposition, the non-unity {\Z\ < 1) renormalization constant Z (which should not 
be confused with Zg appeared in the Hamiltonian) is correctly obtained. This decomposition formula turns out to be 
extremely useful in the following analyses. 

Using these results, we study the nonequilibrium time evolution of interacting bosons in a leaky box (section |rv| ). 
The time evolution is induced if one fixes N at some time (thus the boson state at that time is the NSIB), because in 
most real systems there is a finite probability of exchanging bosons between the box and the environment, hence N 
will fluctuate at later times. We simulate this situation by the following gedanken experiment: At some time t < 
one confines exactly N bosons in a box of volume V, and that at t = a small hole is made in the box, so that a 
small leakage flux J of the bosons is induced. We concentrate on the analysis of the most interesting and nontrivial 
time stage; the early time stage for which Jt <C N, because it is clear that at later times the system approaches the 
equilibrium state. We are interested in the reduced density operator of the bosons in the box; p{t) = Tr^[/5*°*'''(t)], 
where jf°^°'\t) denotes the density operator of the total system, and Tr^ the trace operation over environment's 
degrees of freedom. We successfully evaluate the time evolution of p{t) by a method which is equivalent to solving 
the master equation. Our method gives a physical picture more clearly than the master-equation method. We obtain 
p{t) in a general form in which all the details of the box-environment interaction H^^ have been absorbed in the 
magnitude of the leakage flux J. We show that the time evolution can be described very simply in terms of 60, as 
the evolution from a single NSIB at t < 0, into a classical mixture, with a time dependent distribution, of NSIBs of 
various values of at t > [Eq. ([l07|)]. 

We then discuss the phase (/> as a variable approximately conjugate to the number A^ (section ^). To treat the 
quantum phase properly, we consider the sine and cosine operators, sin <p and cos <j). It is generally very difficult to 
define such operators for interacting many-particle systems. Fortunately, however, in terms of the natural coordinate 
bo we successfully define sin and cos (j) for interacting bosons, using which we can analyze the pha se pr operty in a fully 
quantum- mechanical manner. We define a "coherent state of interacting bosons" (CSIB) [Eq. ( |119| )], which, unlike 
Bogo liubov's ground state |ao,y^')^', exactly ha s the minimum value of the number-phase uncertainty product [Eq. 
(128)]. We also define a new state |^, A, y) [Eq. (133)], which we call the "number-phase squeezed state of interacting 
bosons" (NPIB), which has a smaller phase fluctuation than the CSIB, while keeping the number-phase uncertainty 
product minimum [Eq. (145)]. We point out that p{t) for t > Q can be represented as the phase-randomized mixture 
(PRM) of NPIBs. Among many possible representations of p{t), this representation is particularly convenient for 
analyzing the phase fluctuations and the order parameter. 

We also discuss the action of the measurements (or, their equivalents) of A^ and of (section VI). The forms of 
p{t) after such measurements are discussed. As an example of the phase measurement, we discuss an interference 
experiment of two condensates which are prepared independently in two boxes. It was already established for non- 
interacting bosons that the interference pattern is developed for each experimental run (although the interference 
disappears in the average over many runs) [pT|-p^. Using our formula for interacting bosons, we show that the same 
conclusion is drawn very clearly and naturally for interacting bosons. 

We finally consider the order parameter according to a few typical definitions, as well as their time evolution 
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(section VII). We show that the ofF-diagonal long-range order (ODLRO) does not distinguish NSIB, NPIB and CSIB. 
Hence, the order parameter S defined from ODLRO [Eq. ( [l66|) ], does not distinguish them either. On the other hand, 
the other order parameter VP, defined as the expectation value of the boson operator ijj, has different values among 
these states. In particular, for each element of the PRM of NPIBs, we show that ^ evolves from zero to a finite 
value very quickly: After the leakage of only two or three bosons, each element acquires a full, stable and definite 
(non-fluctuating) value of 'J. 



II. PRELIMINARIES 



A. Phase transitions in finite systems 



We consider the phase transition of interacting bosons confined in a large, but finite box of volume V . Phase 
transitions are usually discussed in systems with an infinite volume (or, the V ^ oo limit is taken at the end of the 
calculation), because infinite degrees of freedom are necessary in the relevant energy scale for strict transitions 

111- In 

such a case, we must select a single physical Hilbert space among many possibilities, which selection corresponds to a 
strict phase transition. However, phase transitions do occur even in systems of finite V in the sense that a single phase 
lasts longer than the time of observation if its linear dimension exceeds the correlation length at the temperature of 
interest [ |l4| . Hence, it is physically interesting and important to explore phase transitions in finite systems. Because 
of the finiteness of V (and the fact that the interaction potential U is well-behaved), the von Neumann's uniqueness 
theorem can be applied. This allows us to develop a theory in a unique Hilbert space. However, since V is large, some 
quantities, which become anomalous in the limit of V — *■ oo due to a strict phase transition, behave quasi anomalously. 
In later sections, we shall identify such a quasi anomalous operator, and discuss how an order parameter is developed. 



B. Effective Hamiltonian 



We start from the standard Hamiltonian for interacting bosons confined in a large, but finite box of volume V: 

Hb= j dVv;t(r) (^-^V^^ V^(r) + iy" <frj d^r'ij\v)i^\v')U{Y - v')^{v')^{v). (1) 

Here, we neglect a confining potential of the box because in the present work we are not interested in its effects such 
as the spatial inhomogeneity of the boson states in the box, and also because we expect that main physics of the 
nonequilibrium evolution of our interest would not be affected by such effects. (Mathematically, our model under the 
periodic boundary condition assumes bosons confined in a three-dimensional torus.) The r dependence of the boson 
field (in the Schrodinger picture) can be expanded in terms of plane waves as 

'/'(r) = -jffaa + ^^k, (2) 



where Op and aj, are called creation and annihilation operators, respectively, of bare bosons. The total number of 
bosons is given by 



N : 



I d\ij\r)i:{r) ^Y.^i^^- (3) 
Jv 



We assume zero temperature, and consider the case where the interaction is weak and repulsive [Eq. (^) below] , and 
where the boson density n is finite (hence, since V is large, N ^ 1): 

n={N)/V>0. (4) 

In such a case, BEG occurs and typical matrix elements of ciq, and N are huge, whereas those of and 
(with k 7^ 0) are small. Taking up to the second order terms in these small quantities, and using the identity. 
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^0^0 + X^k^iio ^-k^-kj ^6 obtain the effective Hamiltonian H in the following form. 
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k#0 \ / 



k^O 



(5) 



Here, e^"-* denotes the free-particle energy, e^^-* = Ti^k^ /2fn, and g is an effective interaction constant defined by 



(6) 



Here, a is the scattering length, and Zg is the first-order "renormalization constant" for the scattering amplitude (hsl 



Z. 



^ ,(o: 



a - 2V ^ M ' 



(7) 



The formal divergence of the sum in Eq. (|^) does not matter because the final results are independent of Zg 
hence the renormalization is successful. We have assumed that 



< na^ < 1, 



(8) 



under which the approximation Hb ~ H is good. We have confirmed by explicit calculations that the term — {.g(l 
Zg) /2V}aQao in Eq. (||) gives only negligible contributions in the following analysis. We thus drop it henceforth: 



N' 



H = 9il + Zg)— + Y.(e^ 



N 



, + q — I a; air 

k^^O 
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2V 



kT^O 



(9) 



Since this H commutes with iV, we can in principle find its eigenstates for which N is exactly fixed. In each subspace 
of fixed N, H is equivalent to 



H{N) ^ + Zg)n^V + ^ e'^aj^Sk 



where 



/ - (0) , 



Note that if we regard aq in H{N) as a classical complex number 

ao -> e*"*V^= "0, 
and al as q;q, we would then obtain the "semiclassical" Hamiltonian iJ"^' as 



^.g(l + Zg)n'V + J2 44«k + ( igne^^^ ^ 4^^^ + h-C 1 
k#o \ k#o / 



(10) 



(11) 



(12) 



(13) 



where we have replaced A^o with N in the last parenthesis because the replacement just gives correction which is of 
higher order in g. This Hamiltonian can be diagonalized exactly (See, e.g., Ref. in which cj) = 0). We shall utilize 
this fact later to find the ground state of H. 



C. Known results for non-interacting bosons 

When 3 = 0, the ground state of free bosons whose number N is fixed, is simply a number state. 
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|iV)^ ^(4)^10), 



/TV! 



(14) 
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where |0) denotes the vacuum of the bare operators; 

OklO) = for all k. (15) 

The energy of \N), 

En = Ne^°^ = 0, (16) 

is degenerate with respect to A^. Hence, any superposition of |A^) is also a ground state. For example, the coherent 
state 

\a) = e-l"l'/2 ^\N) = e-l"l'/'e"'^J|0), (17) 



where a = e*"^-\/]V, is also a ground state that has the same expectation value of N as \N). On the other hand, \a) 
has a finite fluctuation of N, 

(SN')^(iN-{N)r)^\a\'^{N), (18) 
whereas \N) has a definite N. The inverse transformation from |a) to |A^) can be accomplished as 



\N) = I ^\a). (19) 



D. Known results for H 



cl 



Neither \N) nor \a) is an eigenstate when <? > 0. If we can regard oq and a\ as classical numbers ao (= e"^\//Vo) 
and ag, respectively, we can use H'^^ as the Hamiltonian, and its ground state was given by Bogoliubov as ||l|,0,D 



|ao,y ) = exp 



aoa, 



q#0 



where 



cosh \yf 
sinh 



.(0) 



-gn 



2e 
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gn 



Here, is the quasi-particle energy. 



whose dispersion is linear for e'q^ <^ gn; 



The ground state energy is calculated as [|l5| 



^2egie, + 6^°'^ +gn) 
v/.f(.r+2,n). 



1 ^ 128 /na3 



1 



-5n7V + 0(5'-^). 



|0), 



(20) 

(21) 
(22) 
(23) 

(24) 
(25) 

(26) 
(27) 
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The absence of the 0{g^) term in E^^^i means that the large (formally divergent because Zg — > oo) positive energy 
gZguN /2 of H'^^ is canceled by a large negative term arising from the pair correlations. The expectation value and 
variance of N for |ao:y^')^' are evaluated as 

{NY'=N, + Y,{^in\,\y'j^\f, (28) 
(5iVY' = (iV)=' + ^(sinh|y^^'|)^ (29) 



E. Ground state of a fixed number of bosons 

In analogy to Eq. it is possible to construct an approximate ground state of H{N) from |aoiy^')^' as 

\N,y-Y^ r ^l"o,y^y', (30) 

where TV = (A^)"^', Eq. (p8|). We can obtain an explicit form of jA^, y'^')'^' in the form of an infinite series expansion, 
by inserting Eq. ( |2C| ) into the right-hand side (rhs) of Eq. (|30|), expanding the exponential function, and performing 
the (j) integral. However, such an expression is not convenient for the analysis of physical properties. Note that in 
many-particle physics, it is generally difficult to evaluate physical properties even if the wavefunction is known. It is 
therefore essential to find the ground state in the form that is convenient for analyzing physical properties. 

Several formulations were developed for the condensation of interacting bosons with a fixed N. Lifshitz and 
Pitaevskii ]l5[ | developed a formal discussion for the case of fixed N. However, they did not treat ao as an operator, 
hence N was not conserved. For example, 6||m, N) (in their notations) did not have exactly -I- 1 bosons. To treat 
interacting bosons with fixed TV more accurately, one has to include quantum fiuctuations of all modes (including 
k = 0), by treating ao as an operator. Such treatment was developed, for example, by Girardeau and Arnowitt 
Ip^ , Gardiner |^^, and Castin and Dum jl^. A variational form was proposed in Ref. jl^ for the wavefunction 
of the ground state. The variational form takes account of /oMr-particle correlations in an elaborate manner, and 
is normalized exactly. On the other hand, in Refs. jl^ and no explicit form was derived for the ground-state 
wavefunction. [These references are more interested in excited states and the spatially inhomogeneous case, rather 
than the ground-state wavefunction.] We here use a variational form, which is similar to that of Ref. as the 
ground state of H{N); 

\N,y)^e^^(y^^ialr\0). (31) 



Here, G'(y) is the hermite operator defined by 



G(y) ^ ^&Wo E l<l«q«-q + ^-c, (32) 
where y = {yq} are a set of variational parameters, which are taken as 

y, = \yt\- (33) 

Using the well-known formula for arbitrary operators A and B, 

e^Be-^ ^B + [A,B] + ^[A,[A,B]] + ---, (34) 



we find from Eqs. (|l|)-(||) that 



EN,y = {N, y\H\N, y) = ^gnN + o{g^). (35) 



where o{g^) denotes terms which tend to zero as g ^ 0, faster than g^. This demonstrates that the large (formally 
divergent because Zg —^ oo) positive energy gZguN /2 in H{N) is canceled by a large negative term arising from the 



6 



four-particle correlations of the state of Eq. (|3l|). Moreover, we can also show that in the macroscopic limit (V^ — > oo 
while keeping n constant), E^-^y becomes as low as E^^^, of Eq. (^6|). Therefore, the form of Eq. ( ^l|) is a good 

approximation to the ground state. Note that |iV, y) is an eigenstate of N; 

N\N,y) ^N\N,y), (36) 
hence \N,y) is an (approximate) groundstate of H; 

H\N,y) =H{N)\N,y). (37) 



Note also that \N,y) is exactly normalized to unity because e*'^'-^' in Eq. (|3^) is a unitary transformation (although 
it becomes non- unitary in the limit oiV^oo). 



We should make a remark here. In the case of \a,y ) discussed in subsection [ID, E^^^, becomes low enough 



only for the specific choice of the phase of [Eq. (|2l|)] . This phase relation is sometimes called the "phase locking" 
0. From this viewpoint, it is sometimes argued [[7|J19[ that "having a definite phase" and the "phase locking" are 
necessary to achieve a low energy. However, such a statement is rather misleading: In our case, the phase locking 
corresponds to the fact that y^'s are real and positive. On the other hand, our ground state \N, y) has no fluctuation 
in N, hence has no definite phase | pO| because of the number-phase uncertainty relation, 

{SN^){S^^)>l/4. (38) 

Nevertheless, the energy of \N,y) is as low as E^^^i. That is, "having a definite phase" is not necessary to achieve 
the ground-state energy, and thus the term "phase locking" should be taken carefully. 

F. Ground state of A'^ — AA'^ bosons 
The ground state of iV — AN bosons is given by Eqs . (|3l| ) and ( ^2[ ) in which N is all replaced with N — AN. 



However, we are interested in the case where [see section [V 



\AN\ < N. (39) 

In this case, G(y) and y (which are functions of the density of bosons) of iV — AiV bosons are almost identical to 
those of N bosons because {N — AN)/V w N/V = n. Therefore, we can simplify the calculation by using the same 
G(y) and y for all AA^. That is, we take 

|^-AiV,y) :^e'^(y)^= ^ (at)^-^^|0), (40) 

' ' ^{N-AN)r I ^ ^ 

where G{y) and y are those of N bosons. Despite the approximation, this state is exactly normalized and has exactly 
TV - AN bosons; 

N\N - AN,y) = {N - AN)\N - AN,y). (41) 



III. NATURAL COORDINATE 

A. Nonlinear Bogoliubov transformation 

Since we assume that V is finite, e^'^^^^ is a unitary operator (which, however, becomes non-unitary in the limit of 
y — > oo). Utilizing this fact, we define new boson operators 6k by 

6k = e'<^(y)ake"'^(^^ (42) 
This operator satisfies the same commutation relations as dk; 

[^P,^^]=W [bpA] = [blbl] = 0- (43) 
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Note that these relations are exact, in contrast to the operator 6k (k 7^ 0) of Ref. Owing to the exact commutation 
relations, we can define the vacuum of 6k's by 

&k|0,y)=0 for all k. (44) 

From Eqs. (|l|), (||) and (H), we have ^ 

\0,y)^e:^<^y^\0). (45) 

The transformation (^) somewhat resembles the Bogoliubov transformation which diagonalizes i?"^' However, 
in contrast to Bogoliubov's quasi-particles (whose total number differs from TV as an operator), the total number 
operator of 6k's is identical to that of Ok's because [N, G(y)] = 0; 

J2 blb^ = e'^'-y^Ne-'^'^y^ = N. (46) 

k 

This property is very useful in the following analyses. On the other hand, the transformation ( ^ ) is much more 
complicated than the Bogoliubov transformation: The latter is a linear transformation connecting the bare operators 
with quasi-particlc operators, whereas the former is a nonlinear transformation between the bare operators and the 
new boson operators. For example, bo defined by 

bo = e™aQe-"^^y^ (47) 

is a rather complicated, nonlinear function of Ok's and a^^'s (of various k's including k = 0), as can be seen using Eq. 
(^). Such a nonlinear operator bo, however, describes the physics quite simply. Namely, we show that So (and SJ) is 
a "natural coordinate" p,2\ of interacting bosons in the sense that many physical properties can be simply described. 
(It is crucial to find such a coordinate for the analysis of many-particle systems, because generally the knowledge of 
the wavefunction is not sufficient to perform the analysis.) For example, from Eqs. (^), (^) and (Q), we find that 
in terms of 60 the ground state \N,y) is simply a number state; 

\N,y) = ^iblr\0,y). (48) 

In particular, 

bo\N,y)^VN\N^l,y). (49) 

Since 6k 's of k 7^ commute with 60, we also find that 

bk\N, y)^0 for all k 7^0. (50) 

Therefore, in terms of the new boson operators the ground state of the interacting bosons can be simply viewed as a 
single-mode (k — 0) number state. Note, however, that this does not mean that H were bilinear and diagonal with 
respect to 60. In fact, if it were the case then the energy EN,y would be linear in N, in contradiction to Eq. ( |35| ) 
which shows i?Af,y oc N'^ (recall that n = N/V). 

The usefulness of 60 are strongly suggested by Eqs. ( p8| ) and (pO[). We will show in the following discussions that 
this is indeed the case. 



B. Decomposition of 7/) 

Some matrix elements of 60 become anomalously large, among the ground (and excited) states of different A^. For 
example. 



(iV-l,y|6o|A^,y) = ViV = VnF. (51) 

This indicates that in the V — > 00 limit (while keeping the density n finite) 60 does not remain an annihilation operator 
of the physical Hilbert space, signaling that a strict phase transition should occur as V co. Since this anomaly 



8 



should have important effects even for a finite V, it is appropriate to separate 60 from the other terms of i/j. That is, 
we decompose the boson field in a finite system as (c/. Eq. (||)) 

V'--=^&o + V'', (52) 



where Z is a complex renormalization constant. Since we have specified nothing about ^' at this stage, this decom- 
position is always possible and Z is arbitrary. Following Ref. , we define the "wavefunction of the condensate" S 
by 

E = (N -l,y\ij\N,y). (53) 
Since S is independent of r (because both jA^ — l,y) and \N,y) have the translational symmetry), we can take Z as 

E = Z'/^N-l,y\bo\N,y)/Vv. (54) 

That is, from Eq. (||), 

^1/2 = E/V^. (55) 
Then, by taking the matrix element of Eq. ( |5^ ) between |iV — A7V, y) and \N,y), we find 

{N - AN,y\il>'\N,y) =0 for ^AiV (|A7V| < iV). (56) 
We now define two number operators by 

N' = J dVv;'t(r)?A'(r), (57) 
No = N-N'. (58) 
Then, from Eqs. (||), (||), (H), and (||), we find 

{N, y\N\N, y) = "t^jSp + (TV, y\N'\N, y). (59) 

Hence, from Eq. (|58|), 

{N,y\No\N,y) ^V\E\\ (60) 

which may be interpreted as the "number of condensate particles" That is, in agreement with the standard 

result |15|, |E!p is the density of the condensate particles; 

\E\^ ^{No)/V = no, (61) 

where we have denoted the expectation value simply by (■••). We can therefore write S as 



S = V^e^-^. (62) 
Wc thus find the formula for the decomposition of tp as 

i' = e''^\f^bo + ^', (63) 
V nV 

which is extremely useful in the following analysis. Note that we have obtained the finite renormalization; 

\Z\ = no/n < 1. (64) 
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C. Relation to the previous work 



Lifshitz and Pitaevskii p5[ | introduced an operator S that transforms an eigenstate with N bosons into the corre- 
sponding eigenstate with — 1 bosons, without giving an explicit form of S. (They de fined 5 through its matrix 
elements between eigenstates with different values of N. However, as mentioned in section HE , they did not give the 
forms of the eigenstates of fixed iV.) They decomposed ip as (Eq. (26.4) of Ref. [|5|) 



ip = E + ip'. 



(65) 



In the present paper, from Eqs. (47) and (pSf), we obtain the explicit expression for S as 



(66) 



Erom Eqs. ( ^ ) and (^2|), we confirm that 



S|iV,y) -S|iV-l,y), 



(67) 



which was assumed in Ref. ||l5| . The operator S characterizes the condensation by having a finite matrix element [ [15[ . 
In the following we will reveal a striking property of S (or, equivalently, 6o); it is a "natural coordinate" of interacting 
bosons. 

On the other hand, other operators, which also characterize the condensation, were introduced in Refs. ||l0| , [l7| , p^ . 
Girardeau and Arnowitt [0 defined (3o = cloNq^^'^, Gardiner introduced A = ao{No/ N)^^"^ (which is an operator 
form of Eqs. (9) and (10) of Ref. |l|]), and Castin and Dum (l8| introduced a$^^ = / dr^l^{r,t)^p{r,t). These 
operators are totally different from S or bo because complicated many-particle correlations, which are included in S 
and 6o, are not included in /3o, A and a$^^. For example, a$^^ is a linear combination of annihilation operators of 
free bosons, whereas bo is a nonlinear function of both the annihilation and creation operators of free bosons. As a 
result, in contrast to Eqs. ( ^9| ) and (|67|), application of either Po, A or a$^^ to the ground state of N bosons does 
not yield the ground state of iV — 1 bosons; it yields an excited state which is not an eigenstate. Moreover, they are 
not a natural coordinate of interacting bosons in the sense explained in the following. Therefore, we do not use these 
operators, although they would be useful in other problems. 



D. Low- lying excited states 

Excited states of a fixed number of interacting bosons were discussed in Refs. pO| , pT| , | l8[ . In the present formulation, 
we may obtain low-lying excited states by the application to \N,y) of functions of 6]j^'s with k ^ 0. However, since 
we do not need any explicit expressions of the excited states in the following analysis, we do not seek for them in the 
present paper. 



IV. TIME EVOLUTION OF BOSONS IN A LEAKY BOX 



The time evolution of a condensate(s) in an open box(es) was discussed previously for the cases of non-interacting 
and for the case of two-mode interacting bosons in Ref. p^ . In the present paper, using bo, 
we study the case of infinite-mode interacting bosons. 



bosons in Refs. pi- 13 



A. Gedanken experiment 

In most real systems, there is a finite probability of exchanging bosons between the box and the environment. 
Hence, even if one fixes N at some time, N will fluctuate at later times. Namely, the boson state undergoes a 
nonequilibrium time evolution when its number fluctuation is initially suppressed. To simulate this situation, we 
consider the following gedanken experiment [Fig. |l|. 
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J 



t 



FIG. 1. Our gedanken experiment. TV bosons are confined in a closed box for t < 0. At t = a small hole is made in the 
box, so that a small leakage flux J is induced, and the expectation value {N{t)) of the number of bosons in the box decreases 
with time. 

Suppose that bosons are confined in a box which is kept at zero temperature, and that the wall of the box is not 
permeable to the flow of the bosons, i.e., the probability of boson's permeating through the wall within a time scale 
of our interest is negligible. If one measures the number of the bosons at a time t = tp (< 0), and if the box is kept 
closed until i = 0, then the density operator p{t) of the bosons in the box for tp < t < is 



Assume that this box is placed in a large room, which has no bosons initially. Suppose now that at i = he makes 
a small hole(s), or slightly lowers the potential of the wall of the box, so that a small but finite flow J of the bosons 
from the inside to the outside of the box becomes possible for i > 0. We study the time evolution for t > of the 
density operator p{t) of the bosons in the box. 

The expectation value of N will be a decreasing function of t, which we denote {N{t)) (hence (A^(0)) — N). It is 
obvious that as t —^ oo the system approaches the equilibrium state. Therefore, we are most interested in the early 
stage of the time evolution, for which 



Note that if J were not small enough, then the state in the box would evolve into a nonequilibrium excited state. The 
property of such a nonequilibrium state would depend strongly on details of the structures of the box and the hole 
or wall. In the present paper, we are not interested in such structure-sensitive states. Therefore, we assume that J is 
small enough that only transitions between the ground states for different values of N are possible. 



= |7V,y)(iV,y| fortp<<<0. 



(68) 



\N-{N{t))\ <iV. 



(69) 



B. Total Hamiltonian 



Let V denote the volume of the room, which is much larger than the volume V of the box; 



V > K 



(70) 



The total boson field 



^totai^j.^ is defined on V, 



^*°*-'(r) =^(r)+^^(r). 



(71) 
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where ip{r) is localized in the box, and = V'*"*"'— -iA is the boson field of "environment." Then, the total Hamiltonian 
may take the following form: 



where H and are the Hamiltonians of the V' and fields, respectively, when they are closed. Small but finite 
amplitudes of scattering between the i/j and fields are caused by the residual Hamiltonian H^^ . Under our 
assumption (^), the probability of multi-particle collisions during the escape from the box is negligible. Therefore, 
H^^ should take the following form: 

if^^ = A / d3rVi^t(i.)/(r)^(r) + h.c, (73) 



where A is a constant which has the dimension of energy, and /(r) is a dimension-less function which takes values of 
order unity when r is located in the boundary region between the box and the environment, and /(r) = otherwise. 
Although the value of A and the form of / depend on the structures of the box and the hole or walls, our final results 



(e.g., Eq. ( 107)) are independent of such details. 



C. Low- lying states of the total system 

We here list states of the total system which are relevant to the following analysis. Since H^^ is weak, quasi 
eigenstates of the total system are well approximated by the products of eigenstates of the box and of the environment. 
Recall that |iV — {N{t)) | <C iV for the time interval of our interest, and that J is small enough so that only transitions 
between the ground states for different values of N are possible. Therefore, among many possible states of the box 
the relevant states are — AiV, y)'s with |AiV| <C N. On the other hand, there are no bosons in the environment 
at t < 0. That is, the environment is initially in the vacuum, which we denote |0^). Hence, from Eq. (|68|), the initial 
density operator of the total system is 

p*°*'^'(t) = |7V,y)|0^)(0^|(7V,y| {t<0). (74) 

Bosons escape from the box into the environment for t > 0. Since V S> V^, the boson density of the environment is 
kept essentially zero, and BEG does not occur in the environment, for the time period of Eq. (|69|). We can therefore 
take the simple number states \n^,n^r, ■ ■ •)'s of free bosons as eigenstates of the environment, where rik denotes the 
number of bosons in mode k. For example, we shall write |1^) to denote the environment state in which mode k are 
occupied by a single boson whereas the other modes are empty. Therefore, the relevant states of the total system, 
i.e., low-lying quasi eigenstates of can be written as 

\N-AN,y)\n^,n^„---) {\AN\IN), (75) 
where, since conserves the total number of bosons, 

AiV-5]nf (76) 

k 



D. Time evolution in a short time interval At 



We are interested in the reduced density operator of bosons in the box for t > 0; 

p{t) = TrE[p*°*"'(i)], (77) 

where Tr^ is the trace operation over the environment degrees of freedom. The expectation value of any observable 
Q in the box can be evaluated from p(t) as 

{Q{t)) EE Tr*°*^n/5*°*'''(i)Q] 

= Tr[/i(t)g], (78) 
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where Tr denotes the trace operation over the degrees of freedom in the box. Equation ([74|) yields 

m^\N,y){N,y\ 

{Qm = {N,y\Q\N,y). 



(79) 
(80) 



Although we may evaluate p{t) by solving a master equation, we here present a different (but equivalent) method, 
by which the physical meaning can be seen clearly. We begin with noting that a single action of H^^ of Eq. (|7^) can 
only change (A^, N^) by either (— 1,+1) or (+1,-1), and that the latter change is impossible for the initial density 
operator (|74|). Therefore, after a short time interval At which satisfies J At <C 1, the state vector jA^, y)|0^) evolves 
into a state of the following form: 



-iBjv.yAt/ri 



|A^,y)|0^) + ^4')(At)e-'(^"-i-+^i°')^*/^|A^ - l,y)|l^) + 0(A2 



(81) 



where 



(At) 



1 



At 



ifi Jq 

Mk= (lf|(Af-l,y|iJ^^|7V,y)|0^), 

/i = EN,y — -EjV-l,y 

Therefore, the reduced density operator is evaluated as 

p{At) = w{0; At)\N, y) {N, y\ + w{l; At)\N - 1, y) (A^ - 1, y| + OiX^), 

where 

u;(0;At) = l-w(l;At), 
w{l;At)^Y.\^l'\At)\'. 



(82) 

(83) 
(84) 



(85) 

(86) 
(87) 



Here, we have normalized p{At) to order A by Eq. (p6[). We now take At in such a way that 



h/Ec < At < 1/ J, 



where Ec is the energy range of e'"^^ in which |Mkp is finite and approximately constant. Then, since k of the 
environment takes quasi-continuous values, w(l; At) becomes proportional to At: 

w(l;At) = JAt. (89) 

To evaluate J, we calculate |Mkp using Eqs. ( p6|) and (|73| ) as 



A / dV/(r)(l^|(A^-l,y|V'^t(r)s|A^,y)|0^) 



J d\f{v)v^*{v) 



(90) 



where is the mode function of mode k of the environment. Regarding the volume dependence, behaves as 

^ 1/ "v/V, whereas / is localized in the boundary region, whose volume is denoted by t>, of the box and the environment. 
Therefore, 



n VV 



(91) 



On the other hand, from Eq. (pq), an escaping boson has the energy of gn. The density of states of the environment 
at this energy is 



V 

2^ 



m?gn 



(92) 
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Therefore, the leakage flux J is estimated as 



^-iV^^^^y^-^xA^, (93) 
n V h fi 

where numerical factors of order unity have been absorbed in A. We observe that J is reduced by the factor v/V 
1), which means that the escape process is a "surface effect", i.e., it occurs only in the boundary region. On the 
other hand, J is enhanced by the factor N 1). This enhancement is typical to the boson condensation. 

From our assumption of small J, the rhs of Eq. ( |93| ) should be smaller than the critical value Jf-r of the flux above 
which bosons in the box get excited [^. It is seen that this condition is satisfied when v and/or |A| is small enough. 



E. Time evolution for < Jt < iV 



We have found in the previous subsection that the reduced density operator p evolves from the pure state ( [79| ) to 
the mixed state ( |85| ) after a small time interval At (<C 1/^)- Since the latter is a classical mixture of two different 
states, |A^, y) (A^, y| and — 1, y){N — 1, y|, we can separately solve the time evolution for each state. For each state, 
further transitions occur in the subsequent time intervals, (At, 2 At], (2 At, 3 At], • • •. Since V ^ 1^, the recursion time 
of an escaped boson to return to the original position is extremely long (except for rare events whose probability — s- 
as V — > oo.) Therefore, as long as J ^ Jen we can neglect any quantum as well as classical correlations between 
transitions of different time intervals. This allows us to take the no-boson state, |0^), as the initial state of the 
environment for each time interval (^At, (£ + l)At], where ^ = 0, 1, 2, • • •. Hence, for every time interval, we may use 
the same formula (|85|). Furthermore, we can neglect the N dependencies of w and J under our assumption of Eq. 
(H). Therefore, 



-i,y|} 

2,y)(7V-2,y|} 



/5(2At) 

= w{Q- At) {u;(0; At)|7V, y) {N, y| + w{l- M)\N - 1, y) {N - 
+«;(!; At) {u;(0; At)|7V - 1, y){N - 1, y| + w{l- At)|iV - 

= «;(0;2At)|iV,y)(iV,y| 

+«;(!; 2At) |iV-l,y) (TV -l,y| 
+«;(2;2At)|iV-2,y)(iV-2,y|, 

where 

u;(0;2At) = w(0; At)^ 

= (1 - JMf 
w{l] 2At) = w(0; At)w(l; At) + w(l; At)w(0; At) 

= 2(1 - JAt)JAt 
w{2; 2At) = w{l; At)w{l; At) 

= j2(At)2. 

The time evolution in the subsequent times can be calculated in a similar manner. Let 

t = AfAt, 

where M (< N) is a positive integer. We find 

M 



p{t) = ^ w{m;t)\N - m,y){N - m,y| 



(94) 

(95) 
(96) 
(97) 

(98) 
(99) 



where w{m;t) is the binomial distribution; 



Af 



w{m-t) =[ (1 - JAt)*'-™(JAt)^ 



(100) 



We find from Eq. (Q) that w{m] t) is the probability of finding N ~ m bosons in the box at t. From the conservation 
of the total number of bosons, Eq. ([76|), this probability equals the probability that m bosons have escaped from the 
box by the time t. Using Eq. (100), we find 
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M 



{N{t)) = Tr[p{t)N] = ^ w(m; (TV - to) = TV - Ji (101) 

m— 
M 

{N'^it)) = Ti[p{t)N^] = ^ w(m; i)m = J. (102) 

m=0 

Since Ec in Eq. (|8^) is of the order of the atomic energy, we can take At extremely smah such that 



il/ > 1 and M > Jt, 



for a finite t that satisfies Eq. (|69|). In this case, Eq. ( |100| ) can be approximated by the Poisson distribution; 



w{m;t) w K{M, Jt)e- 



Here, K is the normahzation factor. 



1 



K{M, x) 



EX 



m=0 



which approaches unity for aU x as Af — > oo. For large but finite M, we can easily show that 

M+l 



1 



K{M, x) 



1 



v/27r(M + 1) V^'-^ + l 



ex 



(103) 



(104) 



(105) 



(106) 



Therefore, under t he c ondition ( |103| ), we can take K{M, Jt) = 1 to a very good approximation, and we henceforth 
drop K from Eq. ( |104| ). Furthermore, since w{m,t) « for to ^ Jt, we may extend the summation of Eq. ( p9[ ) to 
N. We thus obtain 



N 



pit) « e-^* 

711 — C 

N 

E 



TO 



■|Af-TO,y)(iV-TO,y| 



771 — 

^ ( Jt)^-" 



m=0 



(A^-to)! 



'7^,y)(m,y| 



(107) 



Since this final result is valid even at t = (despite our use of the assumption M ^ 1), it is valid for all t as long as 

< Jt < iV and TV > 1. (108) 

Note that the final result ( |107| ) is quite general because all the details of the box-environment interaction H^^ have 
been absorbed in J. 

The probability P{m, t) of finding to bosons in the box at t is evaluated as 



P{m,t) = w{N - m;t) 



{N - to)! 



(109) 



We call this distribution the "shifted Poisson distribution" , because it is obtained by shifting the center of the Poisson 
distribution. The expectation values and variances are evaluated as 



{N{t)) ^N-Jt 
(N^it)) = Jt 

{SN{t)^) = (SN^it)^) = Jt. 



(110) 
(111) 
(112) 
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V. NUMBER VERSUS PHASE 



A. Cosine and sine operators of interacting many bosons 



Roughly speaking, the conjugate observable of the number is the phase. More precisely, however, physical observ- 
ables are not the phase itself, but the cosine and sine of the phase. Namely, any physical measurement of a phase 
actually measures the cosine or sine of the phase In the case of a single- mode boson (i.e., a harmonic oscillator 
of a single degrees of freedom), it has been discussed that various definitions are possible for the cosine and sine 
operators p^ . This ambiguity does not matter in our case, because we are treating the case where the number of 
bosons is extremely large, whereas differences among different definitions appear only when the number of bosons is 
small. On the other hand, the crucial point in our case is how to select a single "coordinate" (dynamical variable) 
with which the phase is defined, among a huge degrees of freedom. To find such a "proper coordinate" is generally 
very difficult in many-body interacting systems. 

Fortunately, we find that bo is the proper coordinate of interacting bosons, with which we can successfully define 
the cosine and sine operators of interacting many bosons by 



sm0 = 



bp + bl 



2\/blbo + l 
1 



--bo - bl 



2\/blbo + l 
1 



2i\blbo + 1 



(113) 
(114) 



2«V^i^o + 1 



These are the same forms as those of a single harmonic oscillator |25|. In our case, however, there are a huge degrees 
of freedom with mutual interactions. As a result, the Hamiltonian does not take the simple bilinear form with respect 
to 6o, hence the motion of bo + b\ is not that of a harmonic-oscillator coordinate. Nevertheless, many formulas for 
the single mode case are applicable if they are based only on the commutation relations of a boson operator. In 
particular, owing to Eq. (pOj), we can treat any states that have the form of Cm\N — m, y) as if we were treating 
a single-mode problem. 

It will turn out in the following discussions that the above operators give reasonable results for the quantum phase 
of interacting bosons. 



B. Number and phase fluctuations of |A'^, y) 



As we have shown in section [II, the ground state of a fixed number of bosons |-/V, y) can be represented simply 
as a number state if we use a "natural coordinate" Sq. Note that this state has a finite fluctuation of ajao due to 
many-body interactions. Nevertheless, the total number of bosons N has a definite value; 



(7V)Ar.y= (7V,y|7V|7V,y) =iV, 
{5N^)N,y={N,y\5N^\N,y)^Q. 

On the other hand, using the simple representation (^8|), we can easily show that 



(cos(/.>Ar,y = (iV,y|cos(/.|7V,y) = 
(sin(/.)Ar,y = (iV,y|sin0|7V,y) = 0. 



(115) 
(116) 



(117) 
(118) 



Therefore, the ground state of a fixed number of bosons does not have a definite phase ||22l, as expected from the 
number-phase uncertainty relation, E q. (p8| ). It was sometimes argued that although N is definite the fluctuation of 
ajoo might allow for a definite phase ||26| . However, our results (117) and (IIS) show explicitly that the fluctuation 
of djoo does not help to develop a definite phase. Note that this is not due to our special choice of the cosine and sine 
opera tors , beca use t he same conclusion is obtained also when the cosin e and sine operators of do are used instead of 
Eqs. (113) and (114). We will touch on this point again in section VII . 
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C. Coherent state of interacting bosons 



We define a coherent state of interacting bosons (CSIB) by 

|«,y) ^e-HV2^ « l^^y)^ (119) 
— n ym\ 

m—O ^ 

which is labeled by y and a complex number, 

a = e"^VN. (120) 

The inverse transformation is 



\N,y) = / ;r^|«,y). (121) 

Regarding the number and phase fluctuations, we can easily show that |q;, y) has the same properties as a coherent 
state of a single-mode harmonic oscillator |^ . Namely, 

(7V)„,y = (a, y|iV|a, y) = \af = N, (122) 
{SN^)^^y = {a,y\SN^\a,y) = \a\^ = N, (123) 

and, for jap = iV > 1, 



= [l-l/(8|a|2) + ...]sin0, (124) 
{Ssin^ 4>)a,y = (Q;,y|(&in(/))^|Q;,y) 

= (l/4|an(l-sin2 0) + ..., (125) 

and similar results for the cosine operator. It is customary to express the results for the sine and cosine operators 
symbolically as 

(0)a,y « 0, (126) 

{54>^)^.y « l/(4|an = l/(4iV). (127) 

Therefore, |a,y) is the minimum-uncertainty state in the sense that it has the minimum allowable value of the 
number-phase uncertainty product (NPUP) [Eq. (^8|)]; 

{5N^)^,y{54>'')^,y^l/A. (128) 

The magnitude of the number fluctuation is conveniently measured with the "Fano factor" _F, which is defined by 

F = {5N'^)/{N). (129) 

For |a, y), we find 

Therefore, using bo, we have successfully constructed a very special state of interacting bosons, \a,y), whose Fano 
factor is exactly unity, and which has the minimum allowable value of the NPUP. This should be contrasted with 
Bogoliubov's ground state |aoi y'^')'^': for which 



2\cl V ^.rsmh l7,'='n4 

^ > 1, (131) 



and the NPUP is larger than 1/4. The CSIB should not be confused with |ao,y^' 



l\cl 
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D. Number-phase squeezed state of interacting bosons 



We define a new state |C,-/V, y) by 



le, N, y) EE VK{N,m e-l«l ^ )- ^ |m, y) (132) 

.✓ifWm.H./.g_^_(i,.,..|„.,,, (133, 
which is labeled by y and a complex number, 

e = e'1e|. (134) 

We henceforth assume that 

ICP < iV and iV > 1, (135) 

which allows us to set K(N, |^p) = 1 to a very good approximation. 

Th e probability P{m) of finding m bosons for the state A^, y) obeys the shifted Poisson distribution [cf. Eq. 

u-,2 |f|2(^f-'^) 

The number fluctuation and the Fano factor are evaluated as 

{N)^,N.y = N, y|iV|C, N,y)=N- (137) 
{SN^)^,N,y = {C,N,y\SN'\^,N,y) = (138) 



As compared with Eqs. (123) and ( |130| ), we observe that the state |£, TV, y) has a very narrow distribution of the 
boson number. On the other hand, |^, N, y) has a well-defined phase po| when 

1 < \^\^ < N. (140) 

In fact, under this condition we can easily show that 

(sin(/))^,jv,y = (C,iV,y|sin(/)|^,7V,y) 

= [l-l/(8|en + ---]sin</. (141) 
{S sin2 </,)^,jv,y = {^,N,y\ (fein 0)2]^, TV, y) 

= (l/4|Cn(l-sin2 0) + ..., (142) 

and similar results for the cosine operator. As in the case of the CSIB, we may express these results symbolically as 

(0)^,jV,y ~ (143) 
{H')^.N.y « l/im')- (144) 

Therefore, just as \a,y) does, |C, A^, y) has the minimum value of the NPUP; 

{5N^)i^NA^4'^)i,N,y ~ 1/4 (for 1 « « N). (145) 

Since each component of the product satisfies {SN^)^,N,y ^ {SN'^)a,y and (^0^)j,Ar,y ^ {6(j)'^)a,y, |^, A^, y) is obtained 
by "squeezing" |a, y) in the direction of N, while keeping the NPUP minimum, (c/. The conventional squeezed state 
has a larger NPUP pq|.) We thus call |^,iV, y) the "number-phase squeezed state of interacting bosons" (NPIB). 
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E. Phase-randomized mixture of number-phase squeezed states of interacting bosons 



We now take 



^ = e'^^/Jt EE e(t) 



(146) 



That is, = Jt. T hen , inequahties ( |135| ) are satisfied because of our assumption (108). We can show by expHcit 
calculation that Eq. (|107|) can be rewritten as 



P{t) 



^ ^\e''^VTt,N,y){e''^VTt,N,yl 



(147) 



Therefore, the boson state in the box can be viewed either as the shifted Poissonian mixture, Eq. ( |107| ), of NSlBs, or 
as the phase-randomized mixture (PRM), Eq. ( 147 ), of NPIBs. Both representations are simply described in terms 
of bo. In contrast, the same p{t) would be described in an very co mplic ated man ner in terms of bare operators. 

We have thus obtained double pictures (or representations) , Eqs. ( |107| ) and (147), for the same physical state p7| , p8t . 
According to the former picture, the state of the box is one of NSIBs, for which the number of bosons is definite (but 
unknown), whereas the phase is completely indefinite. According to the latter picture, on the other hand, the state 
is one of NPIBs, for which the number of bosons has a finite fluctuation {SN'^) « Jt, whereas the phase is almost 
definite (but unkno wn) {S<f)^) ~l/(4Ji). What allows these double pictures is the superposition principle p^ . 

In addition to Eqs. (107) and (147), there are many other ways to express p{t) as different mixtures. Among 
them, Eq. ( |107| ) is t he form in which each element of the mixture has the smallest value of the number fluctuation, 
whereas in Eq. (147) each element of the mixture has the smallest value of the phase fluctuation. Therefore, the latter 
representation is particularly convenient for discussing physical properties that are related to the phase, as will be 
shown in sections VI B and VII . 



F. Origin of the direction of the time evolution 

As the time evolves, the number of bosons decreases as 

{N{t)) -Jt. 



(148) 



As a result, the energy of the bosons in the box decreases with t. For example, for each element of Eq. (147), 

(e*"*yji, TV, yliJle'"^ VJi, TV, y) < (iV, y\H\N, y) for t > 0. (149) 

However, we note that this energy difference is just a consequence of the difference in (N). Namely, \N,y) and the 
NPIB have the same energy if they have the same value of (N) ; 



J4> 



Jt,N + Jt,y\H\e'^VTt,N + Jt,y) « {N,y\H\N,y). 



(150) 



We can therefore conclude that the direction of the time evolution, from \N,y) to the PRM of |^,iV, y), is not 
determined by an energy difference. Hence, it must be due to difference in the nature of the wavefunctions. The 
study of such a nature, however, needs many pages of analysis, which is beyond the scope of this paper, thus will be 
described elsewhere Eal. 



VI. ACTION OF MEASUREMENT OR ITS EQUIVALENCE 



In the previous section, we have obtained the double pictures, Eqs. (107) and (147). Depending on the physical 
situation, either picture is convenient p8]. To explain this point, we discuss two examples in this section. 
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A. Number measurement 



Supp ose t hat one measures N of the boson system whose density operator is given by Eq. (107), or, equivalently. 



by Eq. (147). In this case, the former expression is convenient. In fact, if the measurement error 

SNerr < y/{SN{t)^) = Vli, (151) 

and if the measurement is of the first kind IBOtl, then the action of the measurement is to narrower the number 



distribution as smaU as SNerr, of the rhs of Eq. ( 107 ). That is, the density operator immediately after the measurement 
is generally given by pi| 



N 



Pj^{t)=J2w{m-N)\m,y){m,y\. (152) 

Here, N is the value of obtained by the measurement, and is a smooth function which has the following properties: 

W{m- N)>0, (153) 
W{m-N)^0 for \m- N\>SNerr, (154) 

N 

J2w{m^ N) = 1. (155) 

The detailed form of W depends on the detailed structures of the measuring apparatus, and thus is of no interest 
here. In an ideal case where SN^rr ^ 0, becomes Kronecker's delta, and 

p^{t)^\N,y){N,y\. (156) 

On the other hand, if the measurement is rather inaccurate in such a way that 

SNerr > \/{SN{t)^) ^ Vlt, (157) 

then almost no change of p is induced by the measurement, if it is of the first kind ]3l|] . 

B. Phase measurement 



Suppose that one measures (j> of the boson system whose density operator is Eq. (107), or, equivalently, Eq. (147) 
In this case, the latter representation is convenient. In fact, if the measurement error 



and if the measurement is performed in such a way that the backaction of the measurement is minimum, then the 
action of the measurement is just to fi nd (o r, get to know) the "true" value of 0, to the accuracy of Siperr, among 



many possibilities in the rhs side of Eq. (147). Therefore, the density operator just after the measurement is generally 
given by [5ll 



P4>(t)^ j ^D{<P- 4>)\e''f'VTt,N,y){e^^VTt,N,y\. (159) 

Here, cj) is the value of (j) obtained by the measurement, and D{(j) — ^) is a smooth function which has the following 
properties: 

D((/)-0)>O, (160) 

D((/)-0)«O for \(j3-4>\>6(j)err, (161) 

r ^D(0-^) = 1. (162) 
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The detailed form of D depends on the detailed structures of the measuring apparatus, and thus is of no interest here. 
On the other hand, if the measurement is very accurate in such a way that 

Sft^err < \/ {S(l>^)cN,y = ^j^^ = (163) 

then p will "collapse" into another state whose phase fluctuation is less than (50^)j_7v_y. However, such an accurate 



measurement is practically difficult when Jt ^ I. Therefore, in most experiments we may take Eq. (ISE) for the 
density operator after the measurement, if the measurement is performed in such a way that the backaction of the 
measurement is minimum. 




FIG. 2. Bosons are confined independently in two boxes. The number of bosons in eadi box is fixed for t < 0. If holes are 
made in the boxes ai t > 0, then the leakage fluxes are induced, which exhibit interference. 

For example, suppose that one prepares bosons independently in two boxes [Fig. H, where the number of bosons in 
each box is fixed for t < 0. At i = a sm all h ole is made in each box, and small fluxes of bosons escape from the 
boxes. Since p of each box evolves as Eq. (|l47l), it is clear that these fluxes can interfere at each experimental run 
(as in the cases of non-interacting bosons [lll[-|13| and two- mode bosons j2^) as the interference between two NPIBs 
of two boxes. From the interference pattern, one can m easu re the relative phase (j) of two condensates. After the 



measurement, p of each box would take the form of eq. (15£). Such an experiment would be possible by modifying 



the experiment of Ref. ||3^ in such a way that J becomes small enough. 

VII. ORDER PARAMETER 

The order parameter of BEG can be defined in various ways. For the ground state |Q!,y'^')'^' of the semiclassical 
Hamiltonian H'^\ different definitions give the same result. However, this is not the case for |iV, y) and |^,iV, y). We 
explore these points in this section. 

A. Off-diagonal long-range order 

We first consider the two-point correlation function defined by 

T(ri,r2) = Tr[/5^t (1.^)^(1.2)]. (164) 
The system is said to possess the off-diagonal long-range order (ODLRO) if [p3|-55| 



lim T(ri,r2)7^0. (165) 

This limiting value cannot be finite without the condensation of a macroscopic number of bosons. (Without the 
condensati on, w e simply has lim|i.i_r9|^oo '""(ri, r2) = for the ground state and for any finite excitations.) In this 



sense, Eq. (165) is a criterion of the condensation. 

If the system possesses the ODLRO, it is customary to define the order parameter S by the asymptotic form of T 

as 

T(ri,r2) ~S*(ri)S(r2). (166) 
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According to this definition, we obtain the same results for all of \N,y), iV, y), and |a,y), where ^ = e"^VJi, 
< iV and a = e^l'^N. Namely, using Eqs. (||) and (|6|), wc find 

lim T(ri,r2) = no, hence S = y^e*''', (167) 

|ri-r2|— >oo 

for all of these states. Therefore, neither the ODLRO nor S is able to distinguish between these states. 



B. Definition as a matrix element 



As an order parameter of the state for which N is exactly fixed, Ref. jlj] uses the "wavefunction of the condensate" 
S, as defined by Eq. (|5^). It is clear that this definition is just a special case of that of the previous section. In fact, 
for |iV,y) Eq. (H) yields 



in agreement with Eq. (167) 



S = (iV-l,y|^|A^,y) = V^e'^, 



(168) 



C. Definition as tlie expectation value of i/) 



Another definition of the order parameter is the expectation value of -0, which we denote by 5"; 

*(r) = (V>(r)). (169) 

According to this definition, the ground state | A^, y) of a fixed number of bosons does not have a finite order parameter; 

* = (A^,y|V-(r)|A^,y) =0. (170) 

This result is rather trivial because "0 alters N exactly by one. On the other hand, it was sometimes conjectured in 
the literature |^6| that the expectation value of the bare operator ao might be finite (A^, y|ao| A^, y) 7^ in the presence 
of many-body interactions because the number of bosons in the bare state of k — fiuctuates due to the many-body 
scatterings. However, this conjecture is wrong because by integrating Eq. (17C) over r we obtain 



(A^,y|ao|Af,y) =0. 



(171) 



That is, although ajoo fluctuates in the state |A^, y) it does not lead to a finite (A ^, y|a o|A^, y). 

In our gedanken experiment, p{t) evolves as Eq. (107), or, equivalently, as Eq. (147). For this mixed ensemble. 



* = Tr[p(i)^] = 0. 



(172) 



This is the average over all elements in the mixed ensemble, and corresponds to the average over many experimental 
runs. On the other hand, ^E* of each elemen t, wh ich co rres ponds to a possible result for a single experimental run, is 
differ ent between the two expressions, Eqs. (107) and (147). That is, for each element of the mixtures, \1/ = for Eq. 
(pOT|) because of Eq. (p7C|), whereas for Eq. Wm 



* = (e'-^VJt, A^,y|V'|e*'^\/ji, A,y) 



(173) 



can be finite. As discussed in sections VP and VE , Eq. (147) is the form in which each element of the mixture 
has the smallest value of the phase fluctuation. This indicates that e ach element of the mixture possesses the most 
definite (non-fiuctuating) value of ^' when we take the representation (147), among many representations of the same 
p{t). We are most interested in this case, because is usually taken as a macroscopic order parameter, which has a 
definite value obeying the Ginzburg-Landau equation. 
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FIG. 3; ^Left scale: {(sini;/))^) of the <?!> = element of Eq. (147). The dashed line represents l/(4Jt). Right scale: 1*1/1 defined 

by Eq. (173). The dotted line denotes l^*] = \/rio- Both {(sini;/)) ) and are plotted against Jt, the number of escaped bosons. 



Figure I plots l^-l of Eq. (|17|) as a function of Jt i— the average number of escaped bosons). We find that 
grows very rapidly, until it attains a constant value 



* = (e''^Vjt,iV,y|V;|e''^ 
This value equals ^' of |a, y) with a = e^'^^/N ; 



Jt,N,y) 



for Jt>2. 



a 



(174) 



(175) 



Note here that |^| of \a,y) is renormalized by the factor \/\Z\ = ^Jn^jn because of the many-body interactions. 
We have also plotted ((sin0)^) of the </> = element of Eq. (147) in Fig. |^. This is a measure of the phase fluctuation 



of the = element. Because of the rotational symmetry with respect to 0, we can regard it as a measure of the 
phase fluctuation of every element. We flnd that ((sin0)^) decreases rapidly as Jt is increased, until it behaves as 



((sin(/.)2) « l/(4Ji), 



(176) 



for Jt ^ 3. Therefore, after the leakage oj only two or three bosons, |e"^vJt, iV, y) acquires the full, stable and definite 
(non-fluctuating) values o/^* and 0, and the gauge symmetry is broken in this sense. One might expect that {SN'^) of 
the order of (N) would be necessary to achieve such stable and </> because {SN"^) = (N) for a CSIB. Our result shows 
that this expectation is wrong, because ^E" and (p already become stable when ~ 2, for which {SN'^) = Jt <^ {N). 
Practically, it seems rather difficult to fix N to such high accuracy that 5N ^ 2. In such a case, 5N would be larger 
than 2 from the beginning, and each element of the mixture has the full and stable values of ^ and (f> from the 
beginning. 

We flnally make a remark on the evolution at later times, whereas we have only considered the early time stage for 
which Jt <ti N. It is clear that the system eventually approaches the equilibrium state. However, a question is; what 
is the state after the early stage, but before the system reaches the equilibrium? It is expected that the state would 
be some coherent state. We can show that this is indeed the case p9|: as t oo, p eventually approaches the PRM 
of |a,y), in which \a\^ = {N{t)) (< N) ||] andy is given by Eqs. and (H) with n = {N{t))/V |§. To 

show this, we must extend the theory of section This is beyond the scope of this paper, and thus will be described 
elsewhere poj . 

The summary of the present paper has been given in section ||. 
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